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Abstract 

Transport or kinetic equations are often derived assuming a quasi- 
particle (on-shell) representation of the spectral function. We inves- 
tigate this assumption using a three-loop approximation of the 2PI 
effective action in real time, without a gradient expansion or on-shell 
approximation. For a scalar field in 1-1-1 dimensions the nonlinear evo- 
lution, including the integration over memory kernels, can be solved 
numerically. We find that a spectral function approximately described 
by a nonzero width emerges dynamically. During the nonequilibrium 
time evolution the Wigner transformed spectral function is slowly 
varying, even in presence of strong qualitative changes in the effec- 
tive particle distribution. These results may be used to make further 
analytical progress towards a quantum Boltzmann equation including 
off-shell effects and a nonzero width. 
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1 Introduction 



Current and upcoming heavy-ion experiments at RHIC and LHC have been 
an important motivation for an extensive theoretical study of kinetic or trans- 
port equations beyond the level of the classical Boltzmann equation. In a 
(nonequilibrium) environment the quasiparticle structure of the theory can 
change and for instance the presence of nuclear matter in heavy-ion collisions 
may affect the mass and width of vector mesons (see e.g. [|l| for a discussion 
of the broadening of the p meson). In order to give a consistent treatment 
these effects have to be incorporated in a kinetic description. Also in other 
areas of physics an extension of the classical Boltzmann equation to a more 
general quantum Boltzmann equation is needed. The theory of baryo- and 
leptogenesis requires the solution of a large set of Boltzmann equations with a 
proper inclusion of CP violating and off-shell effects in an expanding universe 
[@]. In Ref. 1^ a description of the dynamics of Bose-Einstein condensates 
using a quantum Boltzmann equation can be found. 

It is well known [§, H, |, that a consistent setup of kinetic equations de- 
rived from a microscopic quantum field theory needs to incorporate two quite 
different aspects: the statistical information typically described in terms of 
one-particle distribution functions on the one hand, and the kinetic infor- 
mation about the spectrum of the theory, encoded in the spectral function, 
on the other hand. There is a possibility to independently approximate the 
statistical and the spectral aspects of the time evolution equations. In the 
derivation of the classical Boltzmann equation one neglects all nontrivial ki- 
netic aspects of the spectral function, which is approximated by an on-shell 
(5-function. A first-order gradient expansion is then sufficient to derive the 
Boltzmann equation, describing the evolution of the statistical information 
only, from the underlying field theory 0, ^, H, . It turns out that relaxing the 
assumption of on-shell quasi-particles makes the problem considerably more 
complicated. Indeed, an analysis of a self-consistent inclusion of a nonzero 
width has been taken up only relatively recently |0] . Similarly, to go beyond 
leading order in a gradient expansion seems to be rather cumbersome 0. 

Recently, a successful numerical solution of the nonequilibrium time evo- 
lution without the use of a gradient expansion or on-shell approximation has 



been presented for a scalar field theory in 1 + 1 dimensions [IC]. Both the 
early-time behavior and the large-time physics of thermalization can be de- 
scribed remarkably well for quantum fields |T^, as well as for the classical 
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field limit put to the test by Monte Carlo methods and numerical integration 



np. The evolution equations can be obtained from the three-loop expan- 
sion of the closed time path generating functional for two-particle irreducible 
(2P/) Green's functions 0, |1^ and comprise the classical Boltzmann equa- 
tion |, |, 0]. The full solution of these equations allows to directly inves- 
tigate the influence of off-shell effects on the time evolution or the validity of 
a gradient expansion. 

In this paper we focus on the kinetic aspects of the evolution equations 
and compute the nonequilibrium time dependence of the spectral function 
without further approximation. The paper is organized as follows. In the 
following section we briefly review the Schwinger-Keldysh formalism for a 
nonequilibrium scalar field. We derive exact time evolution equations start- 
ing from the Schwinger- Dyson equation for the two-point function and use 
the three-loop 2PI effective action to find a consistent approximation to the 
self energy. In Sec. |^ we discuss the initial conditions used in this paper. 
In order to make contact with heavy-ion physics we consider as initial state 
a so-called "tsunami" : a particle number distribution peaked in momentum 
space. Such an initial state is reminiscent of two colliding high energy wave 
packets [1^. We hasten to say that in our (1 -|- l)-dimensional scalar field 
model direct applicability to phenomology is of course limited. We study in 
this paper in particular the time evolution of the Wigner transformed spec- 
tral function suitable for comparison with standard frameworks for kinetic 
equations. Details of the Wigner transform are given in Sec. |. We present 
our numerical results in Sec. ^. In Sec. |^ we make a connection with equal- 
time methods and moment expansions employed before. Our conclusions are 
summarized in Sec. |^. 



2 Time evolution equations 

Spectral function. We consider a real, scalar quantum field theory with 
field operators obeying the Heisenberg equation of motion (D^ + m^)0(x) + 
A0^(x)/6 = 0. We are interested in the nonequilibrium time evolution of the 
spectral function p given by the expectation value of the commutator of two 
fields 

p(x,i/) = z([0(a;),0(y)]_). (2.1) 
where (■) denotes the average over an initial density matrix. Nonequilibrium 
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Figure 1: Keldysh contour C in the complex-time plane. 

initial value problems can be formulated using the Schwinger-Keldysh or 
closed-time-path formalism For a detailed discussion we refer to Ref. [^]. 
In this formalism the theory is formulated along a contour C in the complex- 
time plane, as shown in Fig. 0. Expectation values are defined along this 
contour. The two-point function is usually written as 

G{x,y) = G^{x,y)Qc{x' - y') + G^{x,y)Qc{y' - (2.2) 

with 0-functions defined along the closed time path C and the Wightman 
functions defined as G^{x,y) = {(j){x)(f){y)) = G^{y,x).^ For a real scalar 
field G'>(x, y) = G'<(x, y). In terms of these, the spectral function is given by 

p{x, y) = i[G>(x, y) - G<(x, y)]. (2.3) 

The spectral function can be obtained from the propagator ( p.2|) directly by 
writing 

Gix, y) = F{x, y) - '-p{x, y) [ec(a;° - y^) - Qc{y' - x^)] . (2.4) 

Here we have separated the propagator using the commutator of two fields, 
p{x,y), and the anti-commutator of two fields, the symmetric or statistical 
two-point function 

F(x,y) = i([0(x),0(y)]+) = ^[G>(x,i/) + G<(x,2/)]. (2.5) 

For a real scalar field F{x,y) and p{x,y) are real functions with symmetry 
properties F{x,y) = F{y,x) and p{x,y) = —p{y,x). We use p and F as the 
basic quantities in this paper. 

^ These two-point functions were also used in Ref. pO|. 
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Evolution equations. The full propagator ( |2.4| ) obeys the Schwinger- 
Dyson equation 

G-\x,y)=G,\x,y)-J:ix,y) (2.6) 

with the proper self energy E and Gq^{x, y) = i px+^ra^] 5p'^{x—y). A differ- 
ential equation for G which may be used to study the time evolution for given 
initial conditions can be obtained from ( P^ ) using d'^'^'^z G~^{x, z)G{z, y) = 
5^^^ {x — y) and one finds 

[a + M^x)] G{x, y)+ij d'^+'z s('^-i°^-')(x, z)G(z, y) = -tS^^'^^x - y) . 

(2.7) 

Here we have separated S in a local part and a nonlocal part, 

j:{x,y) = -iS(i°"^')(x) 6^^\x -y) + S('^°"^i°"^')(x, y) . (2.8) 

The local part S*^^"*^^'^ corresponds to a mass shift and we find it convenient 
to write M'^{x) = + S('°=^i)(a;). 

Since our interest focuses on the nonequilibrium dynamics of p we rewrite 
( |2.7|) such that the spectral function appears exphcitly as one of the dynam- 
ical degrees of freedom. Similar to the above decomposition for G we can 
decompose the self energy 

5.{noniocai)(^^^) = Sp(x,y) - E,{x,y) [Ocix' - y') - Qciv' - x')] (2.9) 

into two real functions, the statistical component of the self energy obey- 
ing T,p{x,y) = T,p{y,x) and the spectral component T,p{x,y) = —T^piii^x). 
These self energies are related to the usual S> < as 

Ep(x,?/) = ii[S>(a;,?/) - S<(a;,?/)], 

T.p{x,y) = -[S>(a;,?/) + S<(x,?/)]. 

We also note that the retarded and advanced propagators and self energies 
are related to the spectral components as Gn^x.y) = 6(a;° — y'^)p{x,y) = 
GAiy, x) and Sjj(x, y) = Q{x^ - |/°)Sp(a;, y) = ^^(y, x). 

It is straightforward to express (|2.7| ) in terms of the real functions intro- 
duced above. The resulting equations are 

Jz^ J dzEp{x,z)p{z,y), (2.10) 
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[U, + M\x)\F{x,y) = - J^dz^ JdzJ:p{x,z)F{z,y) 

+ J^z'^ J dz Sf(x, z)p{z, y). (2.11) 

The form of these equations is exact and they are equivalent to the standard 
Schwinger-Dyson identity (|2.6|) . The equations for p and F are exphcitly 
real and causal. A generic feature is the presence of nonlocal "memory" 
integrals. The time integral in Eq. ( |2.11| ) starts at some initial time taken 
to be z'^ = 0. We note that due to the canonical commutation relations the 
spectral function obeys the equal-time properties 



p{x,y)\xO-. 



0, d^op{x,y)\^o=yO = 5'^{x - y). 



(2.12) 



The time evolution of the spectral function is completely determined by ( p.lO| ) 
and ( gl2D . 

Three-loop 2PI effective action. To solve the Eqs. ( |2.10| ) and ( |2.11| ) one 
has to find suitable approximation schemes for the self energies J^f and Sp 
which respect all symmetries and allow to describe the nonequilibrium early- 
time behavior as well as the large-time approach to equilibrium. Time rever- 
sal symmetry can be easily guaranteed by deriving the evolution equations 
from an effective action by a variational principle. A systematic way to 
achieve the latter is by using the loop expansion of the 2PI effective action 
1^, formulated along the closed time path [|15], |[. Recently, it has been 



shown that the equations in the three-loop approximation describe both the 
early-time behavior and the physics of thermalization |TD, TT]. The self en- 
ergies in this approximation |10|, here expressed in terms of p and F, are 
for a Gaussian initial density matrix given by 



^pix,z) 



T,f(x,z) 



-yP(2;,2;) 
--F{x,z) 

D 



F^{x,z) - ^p\x,z) 



(2.13) 
(2.14) 



(cf. Ref. [|T0[ for the extension to quartic initial density matrices). The effec- 
tive mass term including the local part of the self energy reads 



M'^{x) 



m 



-F{x,x). 



(2.15) 
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Eqs. ( p.lO| , |2.11| ) with ( p.l3| - |2rT5| ) form a closed set for p and F that is energy 
conserving and time- reversal invariant. We emphasize that the propagators 
in the "tadpole" diagram contribution in Eq. ( [^.15D and in the "setting- 
sun" contributions (|2.13| , |2.14|) are full ones. The form of the self energies 
presented here (with free instead of full propagators) has also been found in 
a perturbative calculation in thermal equilibrium We note from Eqs. 



(|2.10| , |2.11| ) that in local approximations, such as the Hartree approximation, 
derivable from the two-loop 2PI effective action, and leading-order large N 
schemes, the evolution of the spectral function p does not couple back to the 
evolution of the symmetric two-point function F. 

Thermal equilibrium. It is instructive to consider the introduced quan- 
tities in thermal equilibrium. The two-point function G in equilibrium is 



described by the same equation ( |2.7| ) if the closed time path is replaced by 
an imaginary path C = [0, —i(3], with /3 the inverse temperature [Q. Since the 
equilibrium two-point functions depend only on the relative coordinates it is 



convenient to consider the Fourier transforms p' 



(eq)| 



UJ 



p), etc. From the peri- 



odicity ("KMS") condition for G{x,y) and S(°°°i°^'^i)(a;, |/) in imaginary time 



one finds the generic equilibrium relations relating for example G^^'^\uj,p) 
and gJ^^^cu,?) g 



G 



UJ,p) 



/3a; ^(eq) 



(2.16) 



For the spectral and statistical components we employ here this translates 
into the following equilibrium relations [16 

F'~^^\uj,p) = 



4"^(^,P) 



(2.17) 
(2.18) 



with nB{uj) 



v-1 



While the spectral function p*^'^^-' encodes the infor- 



mation about the spectrum, one observes that the symmetric function F^^'^ 
encodes the statistical aspects in terms of the particle distribution function 



n-Q. For a vanishing u;-dependence the function r(*^'^)(ci;, p) 



p)/2a; 



plays the role of a decay rate for one-particle excited states with momen- 
tum p. 
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3 Initial conditions 



We consider a spatially homogeneous system such that, after spatial Fourier 
transformation, the two-point functions depend only on p and we consider the 
modes p(a:°, y^; p). As indicated above, the initial conditions for ?/°; p) 
are fixed by the equal-time relations (|2.12|) . As initial conditions for the 
symmetric two-point function ?/°; p) we consider a "tsunami" -like situ- 



ation |jT2[ in which a spatially homogeneous collection of particles move with 
approximately the same momentum peaked around pts and — Pts- Explicitly, 
we take 



1 



Fo{x ,y ;p)|^o=j^o=o = — 



0, 



d^odyoFo{x ,y ;p)|^o=yO=o 
9^.oFo(x°,/;p)|^o=yO=o 
with ujp = a/p^ + and the initial particle number 

1 



^o(p) + 



(3.1) 



A/'exp 



2^2 



IPI 



IPtsI)^ 



(3.2) 



Here a controls the width of the initial distribution and A/" is a normalization 
constant. The initial condition is clearly far from thermal equilibrium and 
reminiscent of two colliding high energy wave packets IH 



To study the nonequilibrium evolution of the spectral function p{x , y , p) 



we solve the coupled Eqs. ( p.lO| - |2?TT| ) with ( p.l3| - pT5| ) numerically along the 
lines of ref. [0. We use a standard discretization on a lattice in space 
and time for a finite spatial volume and spatially periodic boundary condi- 
tions. To remove finite size effects we increase the volume until the results 
become stable. Because the evaluation of the memory integrals is rather 
time- and memory- consuming the numerical results presented below are for 
a (1 + 1) -dimensional system. The lattice acts as a cutoff and regulates the 
ultraviolet divergences, which have to be renormalized. In 1 + 1 dimensions 
this is straightforward since only the one-loop contribution is logarithmically 
divergent. We absorb this in a bare mass parameter n with the replacement 



m 



m 



— Sm? in Eq. ( |2.15D . The counterterm 5m? cancels the diver- 



gent vacuum contribution coming from the one-loop graph. The finite part 
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of 5m? is fixed by requiring that the renormahzed one-loop mass parameter 
in vacuum (no(p) = 0) equals m and we express all dimensionful scales in 
units of m. For the plots we have used a space lattice with spacing ma = 0.3 
and a time lattice ao/a = 0.25. The system size is mL = 24. 



4 Wigner transformation 

Derivations of transport or kinetic equations are typically performed by con- 
sidering correlation functions in the so-called Wigner representation. Two- 
point functions are written in terms of the center-of-mass coordinate X = 
{x + y)/2 and the relative coordinate x — y, with respect to which a Fourier 
transformation to momentum space is performed (Wigner transformation). 
Derivations based on a gradient expansion in dx assume a slow variation 
of the two-point functions with X. In particular, the classical Boltzmann 
equation can be obtained from (|2.10| - pA5|) using a first-order gradient ex- 
pansion and an on-shell form of the spectral function with zero width [^]. 
The full equations (|2.10| - prT5|) can therefore be considered as a "quantum 
Boltzmann equation" including off-shell effects and resumming an infinite 
order of derivatives. 

To analyze the evolution of the spectral function (not to solve the dy- 
namics) we perform a Wigner transformation and write 

ip(X°;cu,p) = / dte'^^'piX^ + t/2,X'^ -t/2;p). (4.1) 

J-2X0 

The i is introduced such that p(X°;c<j,p) is real. Because the spectral 
function is antisymmetric, p{X^;uj,p) = — p{X^ ; —u , p) , we will present 
the positive-frequency part only. Furthermore, as a consequence of the 
equal-time relation ( p. 12 ) the Wigner transform obeys the sum rule / 



J duj/{2TT) uj p{X^;uj,p) = 1. In our simulations an evaluation of the sum 
rule typically results in 1.000 < I < 1.003, indicating the accuracy in the 
numerical evaluation of the Wigner transform. 

Since we consider an initial- value problem with > 0, the time 

integral over t = x° — is bounded by ±2X°. This leads to a so-called 
"finite-time effect" in the Wigner transformed quantities. To understand the 
qualitative aspects of a finite-time interval it is instructive to consider the free 
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spectral function po{x^,y^; p) = t^p^ sin[c<jp(x° — y^)] with cOp = a/p^ + m'^- 
After performing a Wigner transformation, we find 

sin[(^-a;p)2X°] sin[(^ + ^p)2X°] 
PolA ,u;,pj- — ■ — . [4:.Z) 

UJp{uJ-UJp) UJp{uJ + UJp) 

For finite X° this spectral function shows a rapidly oscillating behavior, while 
its envelope is peaked at = ±Ljp. In the limit X° ^ oo it reduces to 
po{uj, p) = 27rsgn(u;)5(u;^ - u;^). 

For an interacting system, however, effective damping leads to a finite 
correlation time r, such that finite-time effects are expected to vanish when 
X° ^ T. This can be investigated analytically by assuming that the spectral 
function is strictly exponentially damped and can be written as p{x^, y^; p) = 
g-7p|a; -y \ sm[Ep(x'^ — ?/°)]. The effective mass Ep and rate 7p are al- 
lowed to depend on the time X°. The corresponding Wigner transform 
reads p{X^;u!,p) = pbw(-^°; p) + ^p{^^',^:P), where pew denotes the 
Breit- Wigner function 

PBw(X°; u, p) = (4.3) 



with a width Tp{X^) = 27p(X°). The additional contribution 5p(X°;u;,p) 
vanishes exponentially as exp(— FpX"), such that the finite-time effect is 
indeed absent for large enough (but finite) We want to emphasize that 
this finite-time effect is a general consequence of an initial-time problem 
formulated in terms of Wigner transformed quantities and should be taken 
into account when used to describe the early-time behaviour (FpX° < 1) in 
an experimental setup. 

Finally, we note that the positivity condition, sgn(u;)p(X°; c<j, p) > 0, can 
only be shown to hold in the special case that the initial density matrix 
commutes with the full Hamiltonian, such as in thermal equilibrium. In 
this case the system is of course stationary and independent of X^. As a 
consequence the interpretation of the nonequilibrium spectral function as 
the density of states should be taken with care. 



5 Nonequilibrium evolution 

In Fig. ^ we show the time evolution of the spectral function p(X°+t/2, X^ — 
t/2;p) as a function of the relative time t for two different momentum modes 
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Figure 2: Absolute value of p{X^ /2 + t, /2 — t;p) as a function of t at fixed 
mX^ = 35.1 for spatial momentum p = 27ck/L with A; = (full), 4 (dotted). 
The coupling is A/m^ = 4. 



at given as obtained from the numerical solution of Eqs. ( p.lO| - |2rT5| ). 
From the logarithmic plot one observes an effective damping of the oscil- 
lating spectral function while the oscillations never damp out completely. 
The decrease of the maximum amplitude very quickly approaches an ap- 
proximately exponential behavior characteristic for the decay of one-particle 
excited states with given momentum. One observes that the zero momentum 
mode is stronger damped than the higher momentum mode.0 

The effective damping should correspond to a nonzero width of the Wigner 
transformed spectral function. In Fig. |^ we display the Wigner transforms 
p{X^;uj,p) for the zero momentum mode (left peak) and the higher mo- 
mentum mode (right peak) for the same X°-time as in Fig. |^. One clearly 
observes that the interacting theory has a continuous spectrum described by 
a peaked spectral function with a nonzero width. The inset shows a blow-up 
of the zero mode around the three-particle threshold 3Ep/m = 4.38. The 
expected bump in the spectral function is small but visible. We stress that 



^An analysis of damping times as a function of the coupling A can be found in Ref. 
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Figure 3: Wigner transforms p{X^]UJ,p) of the real-time spectral functions 
of Fig. 1^ as a function of uj at mX° = 35.1 for k = (left peak) and 
k = 4 (right peak). Also shown are fits to a Breit- Wigner function (dotted) 
with {Ep/m,Tp/m) = (1.46,0.37) {k = 0) and (1.82,0.30) {k = 4). The 
inset shows a blow-up of the zero mode around the three-particle threshold 
3Ep/m = 4.38. The expected bump from off-shell decay is small but visible. 

this bump in the spectral function is kinematically forbidden for the on-shell 
limit and arises purely from off-shell decay. In Fig. |^ we also present fits 
to a Breit- Wigner spectral function. While the position of the peak can be 
fitted easily, the overall shape and width are only qualitatively captured. In 
particular, the slope of p(X°;a;,p) for small u is quantitatively different. We 
also see that the Breit- Wigner fits give a narrower spectral function (smaller 
width) and therefore would predict a slower exponential relaxation in real 
time. 

We now turn to the dependence of p(X°;a;,p) on X°. In Fig. ^jwe show 
the spectral function at zero momentum p{X^;uj,p = 0) for three different 
values of X^. We see that the peak slowly shifts to smaller values. In the inset 
to Fig. ||this behavior is quantified by plotting the X°-dependence of Ep{X^) 
obtained from a fit to the Breit- Wigner function ([4.3|) . For comparison, the 
solid line in the same figure shows the time-dependent one-particle energy 
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Figure 4: Time dependence of the spectral function p(X°; tu, p) at three values 
of X° for p = 0. The peak slowly shifts towards lower values. The inset shows 
the X°-dependence of Ep{X^) obtained from a fit to a Breit-Wigner function 
(circles) as well as the one-particle energy ep{X^) defined via equal-time two- 
point functions, see Eq. (|5.2| ) below (line). 

ep(X°) as defined via equal-time two-point functions in Eq. ( ^.2|) below. We 
see a good agreement of both definitions of the one-particle energy for large 
enough times. We come back to this below. 

A qualitative way to characterize the X°-dependence of the shape of the 
spectral function is to assume that it is approximately described by a Breit- 
Wigner function. The time dependence of TplX^) can then be followed by 
performing fits for different values of X^. The resulting rp(X°)/2 is shown 
in Fig. 1^ for two momentum modes. After a variation for early times, the 
X°-dependence becomes weaker for larger times. As indicated before (see 
Fig. the Breit-Wigner fits result in a smaller width than what is expected 
by comparing with the actual spectral function. A simple way to obtain a 
quantitative description is by extracting a relaxation rate jp directly in real 
time, i.e. from a fit of the spectral function p{X^ + t/2,X^ — t/2]p) to an 
exponential form ~ exp[— 7p(X°)t] (see Fig|^). The resulting rates are also 
presented in Fig. ^. We observe that the relaxation rate for the zero mode 
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Figure 5: Time dependence of the relaxation rate 7p(X°), determined directly 
in real time from p(X°+t/2,X°-t/2;p) forp = 27rfc/L, k = (full) and A; = 4 
(dashed). Also shown is rp{X^)/2 from a fit to a Breit-Wigner function, for 
k = (squares) and k = 4 (triangles). 



increases first, corresponding to a broadening of the Wigner transformed 
spectral function, before it levels off at later times. Despite the quantitative 
differences we see a similar qualitative X°-dependence from the Breit-Wigner 
fit for large enough times. 

Particle distribution function. The relatively weak X°-dependence of the 
kinetic aspects of the evolution encoded in the spectral function has to be 
confronted with the time dependence of the statistical aspects encoded in F. 
The particle distribution is expected to exhibit strong qualitative changes, 
since we start with an initial condition ( p.2|) far away from thermal equi- 
librium and the current approximation is known to describe the physics of 
thermalization at large enough times 
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The evolution of the statistical information can be extracted from the 
equal-time correlation function F{t,t]p). If the theory exhibits a quasi- 
particle structure it is possible to define a quasi-particle distribution function 
from a combination of equal-time two-point functions [1^, |1^ . For this pur- 
pose we introduce the two-point function K{x,y) = |([7r(x), 7r(?/)]_|_), where 
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Figure 6: Effective particle number n{ep,X^) as a function of Sp at four 
values of X^. As time goes on, the original "tsunami" vanishes and the 
infrared modes become more excited. 

7r(x) = dxO(f){x) is the canonical momentum. In the spatially homogeneous 
situation we consider here, K{x^,y^;p) = dxodyoF{x^,y^;p). We define the 
distribution function n(£p(X°),X°) at time as 

n{ep, + ^ = [K{X', X°; p)F(X°, X'; p)] . (5.1) 

For the free theory the definition of n{ep,X^) corresponds to the standard 
one-particle distribution. The distribution function depends on a time-dep- 
endent one-particle energy, defined as 

SpiX') = [K{X', p)/FiX', X°; p)] . (5.2) 

The applicability of this definition to describe one-particle energies for the 
present nonequilibrium situation finds an impressive demonstration by com- 
paring it with the position of the peak in the spectral function as discussed 
above (cf. Fig. 

In Fig. ^ we show the results for the effective particle number n{ep, X°) as 
a function of Ep at four values of The initially high occupation number in 
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a small momentum range rearranges as time goes on. The original "tsunami" 
vanishes and the infrared modes become more excited. The distribution is 
not yet thermal. Thermalization takes place on longer time scales which 
is shown for this model in Ref. |jlO[- The strong qualitative changes of the 



statistical aspects encoded in the particle distribution function have to be 
confronted with the comparably moderate changes of the spectral function 
shown in Fig. ^ and ^. 

6 Equal-time spectral function and moments 

In order to make a connection with methods employed before, we discuss 
some equal-time properties of the spectral function. As has been discussed 
earlier, the spectral function is the commutator of fields and obeys elemen- 
tary equal-time properties p.l2| ). There is a further hierarchy of equal-time 



relations that follow from the exact evolution equation ( p.lO|) . The identity 

aXx°,i/°;p)Uo=,o = (6.1) 

can also be traced back directly to the operator equation of motion. The 
following three equal-time derivatives are less trivial and read 



9Xx°,i/°;p)Uo=,o = -[p^ + M\x) 
9^op(x°,y°;p)Uo=,o = -2d,oM\x), (6.2) 

-a^oEp(x°,?/°;p)|^o=yO. 

These expressions are exact. One observes that the self energy Ep, necessary 
to describe scattering, enters for the first time only in the fifth derivative. 
This suggests that in formalisms which are strictly based on equal-time corre- 
lation functions |^, the kinetic aspects of scattering encoded in the spectral 



function may only be captured in very sophisticated truncations. 

The above relations can be immediately recast in terms of moments of 
the spectral function in the Wigner presentation. We define the n'th moment 
as 

/oo J 
-^^>(XO;^,p). (6.3) 
-oo 
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Due to the antisymmetry of p{X^; u, p) all even moments A2n{X^; p) vanish. 
It is straightforward to relate the odd moments directly to the equal-time 
derivatives given above, using d^o = ^dxo + df. The first moment is just 
the sum rule, y4i(X°;p) = 1 and Eq. (|6.1|) provides a consistency check, 
dloPix'^ ,y^',p)\xO=yO = dxoAi{X^;p) = 0. For the higher derivatives one 
finds 

9,%p(x°,i/°;p)Uo=,o = -2dxoAs{X^-p), (6.4) 
d'Mx\y';p)\.o=yO = -^9^oA3(X°;p)+A5(X°;p). 

The combination of ( |6.4|) and (|6.2|) gives information on the time dependence 
of the lowest moments. Therefore, also in approximation schemes that are 
based on an expansion in moments, the lowest quantity that shows a direct 
dependence on the self energy Sp appears only at relatively high (fifth) order. 
We conclude that both in equal-time formalisms and in approximations based 
on moments a proper inclusion of scattering effects might be difficult to 
achieve. 

7 Conclusion 

The presence of a nonequilibrium environment affects the quasiparticle struc- 
ture of the theory, encoded in the spectral function. As a relevant physical 
application one may think of high energy heavy-ion collisions where spec- 
tral properties of the vector mesons are affected by the presence of nuclear 
matter. 

In this paper we have concentrated on a self-consistent determination 
of the spectral function out of equilibrium, using a scalar field theory in 
1 + 1 dimensions as our (toy) model. We have calculated the nonequilibrium 
time evolution of the spectral function by solving the equations of motion 
obtained from the three-loop 2PI effective action in real time. To com- 
pare the dynamics of the spectral function and the statistical aspects of the 
theory we also compute the evolution of an effective particle number distri- 
bution. Scattering and off-shell contributions are included in the dynamical 
evolution, in contrast to well-studied Hartree and leading-order large ap- 
proximations. Even for moderate couplings and initial conditions far away 
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from equilibrium, we observe a rather weak dependence of the Wigner trans- 
formed spectral function on This fact provides a necessary condition for 
a successful gradient expansion typically employed in the derivation of ki- 
netic equations. After variations at early times we observe that the spectral 
function approaches a rather stable shape characterized by a nonvanishing 
width. A consistent "quantum Boltzmann equation" certainly has to relax 
the zero-width assumption employed in the derivation of the standard Boltz- 
mann equation. We have argued that it might be notoriously difficult to 
derive the kinetic aspects of scattering encoded in the spectral function from 
equal-time formalisms or approximation schemes based on an expansion in 
moments. We emphasize that the employed framework based on the 2PI 
effective action allows for a consistent description of both the early and the 
late time behavior. Since a change in the position of the peak and width 
of the spectral function can be described consistently in this framework, it 
would be interesting to extend the study presented here to more realistic 
models for heavy-ion collisions. 

Finally the analysis performed in this paper was largely motivated by the 
physics of relativistic heavy-ion collisions. We would like to stress, however, 
that the approach using the 2PI effective action is general and that possible 
applications of the formalism range from the early universe to condensed 
matter systems. 
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